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1.  This  report  is  the  Final  Report  on  the  contract.  It  covers 
research  done  on  radiative  transmission  lines  during  the  eleven- 
month  period  from  1 Jan  to  30  Nov  1976.  Radiative  transmission 
lines  are  applicable  to  several  radar  and  communication  functions. 

The  particular  electromagnetic  guiding  structure  analyzed  in  this 
report  is  a periodically  slotted  dielectric  coaxial  line  above, 
parallel  to,  or  buriet^in,  a lossy  ground.  The  ground  modifies  the 
propagation  characteristics  of  the  coaxial  line,  reducing  the  range 
at  which  radar  targets  can  be  detected.  The  objective  of  this 
research  is  to  investigate  theoretically  the  effects  of  a lossy 
ground  on  the  transmission  characteristics  of  a radiative  transmission 
line.  A second  problem  considered  is  the  radar  scattering  by 
obstacles,  such  as  vehicles  or  personnel,  in  the  guided  wave  fields. 
The  contractor  carried  out  a rigorous  mathematical  analysis  of  the 
periodically  slott|d  transmission  line  which  relates  its  radiative 
properties  to  its  physical  configuration.  These  relations  can  be 
used  as  the  basis  for  determining  the  effects  of  ground  conductivity 
on  the  performance  characteristics  of  radiative  lines  for  guided 

wave  radar  applications. 

2.  The  above  work  is  of  value  since  it  provides  basic  knowledge  which 
makes  possible  improved  and  optimized  electromagnetic  devices  for 
USAF  security  systems,  radar,  and  communications  systems. 

WALTER  ROT MAN 
Project  Engineer 
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Abstract 


After  formulation  of  an  integral  equation  and  of  the  dispersion 
relation  for  a periodically  slotted  coaxial  cable  above  and  parallel  to 
the  air-lossy  ground  planar  interface,  the  report  addresses  the 
evaluation  of  the  exterior  dyadic  Green's  Function  in  a unit  cell.  The 
latter  constitutes  the  heretofore  not  available  in  the  literature  kernel  of 
the  integral  equation.  The  problem  is  reduced  to  an  evaluation  of  the  field  of 
electric  and  magnetic  multipole  phased  line  sources  above  a lossy  ground.  The 
analysis  of  the  various  relevant  integrals  is  carried  out  either  in  terms 
of  an  asymptotic  expansion  or  rigorously  by  function  theoretical  methods. 
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The  influence  of  lossy  ground  on  the  propagation  characteristics 
of  periodically  slotted  coaxial  cable  is  of  interest  in  guided  wave  radar, 
or  in  intrusion  prevention. 

This  study  presents  an  approach  to  a systematic  development  of 
the  formalism  necessary  for  the  derivation  of  the  dispersion  relation  of 
a periodically  slotted  coaxial  cable  above  and  parallel  to  a flat  lossy  ground- 
air  interface.  Continuity  of  tangential  fields  across  a slot  in  a unit  cell  yields 
an  integral  equation  over  the  slot.  The  kernel  of  the  integral  equation  is  the  sum 
of  the  exterior  and  the  interior  unit  cell  dyadic  Green's  functions,  for  the 
unperforated  coaxial  cable.  This  integral  equation  may  be  solved,  e.  g.  , 
via  Galerkin's  procedure  which  yields  a set  of  homogeneous  linear  equations. 
Setting  the  determinant  equal  to  zero  yields  the  dispersion  relation  for  the 
slotted  cable.  The  exterior  dyadic  Green's  Function  is  separated  into  that 
of  the  cable  in  free  space  plus  a correction  which  contains  the  effect  of  the 
ground.  As  show  n in  Appendix  A (see  also  L 1 ])  , the  exterior  dyadic  Green's 
function  for  the  slotted  cable  in  free  space  may  be  derived  from  two  axial 
Hertz  potentials,  which  are  expanded  in  terms  of  axial  and  angular  harmonics. 
Each  combination  of  a spatial  and  an  angular  harmonic  is  identical  in  form 
with  a Hertz  potential  ofanelectric  or  magnetic  current  phased  line  source 
located  on  the  axis  of  the  coaxial  cable.  Thus,  if  one  neglects  multiple  in- 
teractions, the  effect  of  the  ground  is  obtained  as  a superposition  of  back 
scattered  fields  excited  by  the  various  phased  line  sources.  Hence,  the 
evaluation  of  the  exterior  dyadic  Green's  function  is  reduced  to  the  deter- 
mination of  fields  excited  by  an  electric  or  a magnetic  phased  line  source 
above  and  parallel  to  a lossy  ground.  The  case  of  a phased  electric  line 
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source  was  treated  in  [2].  Reference  [2],  however  .limits  its  formulae 
to  the  case  when  the  dielectric  constant  e of  the  ground  is  sufficiently  large, 
so  that  the  propagation  constant  of  the  ground  may  be  replaced  by  kr/e~  , 

This  procedure  restricts  the  distance  of  the  line  source  from  the  interface 
to  be  sufficiently  large,  so  that  the  higher  order  terms  arising  from  the 
expansion  of  the  ground  propagation  constant  may  be  neglected.  However, 
in  this  case,  the  effect  of  ground  would,  in  general  , be  a small  perturbation 
and  the  propagation  characteristics  of  the  cable  would  be  practically 
independent  on  the  presence  of  ground. 

To  amend  the  situation  for  the  TE  modes  this  report  carries  out  the 
calculations  to  higher  order  in  terms  of  an  asymptotic  multipole  series. 

For  the  TM  case  the  procedure  of  Reference  12]  modifies  the  location  of 
Sommerfeld  poles,  which  in  certain  cases,  particularly  for  not  very  large 
values  of  |s|,  may  give  rise  to  discrepancies.  Furthermore,  the  validity 
of  the  analytic  continuation  based  on  values  of  the  integral  at  four  different 
points  (Appendix  A of  L2])  is  questionable.  On  the  other  hand  the  method  of 
evaluation  given  in  this  report  has  a rigorous  basis;  no  approximation  in  the 
location  of  the  poles  is  made  here  and,  a rigorous  result  is  obtained  via  a 
contour  deformation,  which  explicitly  separates  the  pole  contribution  from 
the  quasistatic  term  plus  a rapidly  convergent  correction. 

The  additional  advantage  of  the  separation  of  the  pole  contribution  is 
that  it  permits  a simple  continuous  tracking  of  this  contribution  in  a numerical 
solution  of  the  dispersion  relation,  when  the  propagation  constant  0 of  the  cable 
becomes  complex.  In  contrast,  the  function  that  appears  in  L2]  may  become 
discontinuous,  a feature  not  admissible  in  a dispersion  relation. 
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The  organization  of  this  report  is  as  follows: 

Section  Z discusses  the  integral  equation  for  a slotted  coaxial 
cable,  and  the  resulting  form  of  the  dispersion  relation  which  is  applii  able 
for  both  the  slow,  and  the  fast  wave  regions. 

Section  3 addresses  the  derivation  for  the  exterior  dyadic  Green's 
function  which  constitutes  the  key  ingredient  in  the  kernel  of  the  integral 
equation.  The  problem  is  reduced  to  the  evaluation  of  certain  canonical 
integrals. 

Section  4 presents  an  asymptotic  evaluation  of  the  TE  contribution. 

Section  5 addresses  the  evaluation  of  the  TM  contributions  in  terms 
of  Somtnerfeld  wave  residue  an  explicit  quasistatic  part  and  two  rapidly 
convergent  integrals.  The  explicit  residue  contribution  permits  a simple 
analytic  continuation  of  the  functions  appearing  in  the  dispersion  relation. 
Appendix  A,  after  summary  of  radial  line  formalism  appropriate  to  unit  cell 
in  a cylindrical  structure,  derives  the  exterior,  unit  cell,  dyadic  Green's 


function  for  the  coaxial  cable  in  free  space.  Appendix  B shows  that  the 
electromagnetic  fields  due  to  phased  line  source  parallel  to  an  interface 
between  two  homogeneous  dielectric  half  spaces  is  derived  from  two  axial 
Hertz  potentials.  The  expressions  for  the  potentials  are  derived.  Appendix  C 
discusses  the  location  of  Sommerfeld  poles  arising  in  the  TM  contribution  to 
the  field  excited  by  line  sources  in  the  presence  of  les  sy  ground. 


2.  Analysis  of  the  Dispersion  Relation  for  a Periodically  Slotted 
Coaxial  Cable  AboveLossy  Ground. 

Figure  (1)  shows  the  geometry  of  a dielectric -loaded, periodically 
slotted  coaxial  cable  with  a zero  thickness  shield.  A dielectric  jacket  is 
not  included,  but  its  presence  could  be  easily  incorporated. 

Figure  (2)  shows  such  cable  above,  and  parallel  to,  the 
planar  interface  between  air  and  absssy  ground.  A unit  cell  of  the  structure  is 
contained  between  two  parallel  planes  perpendicular  to  the  cable  axis  and  a 
distance  d apart.  The  dispersion  relation  for  the  periodically  slotted  cable 
is  obtained  with  the  help  of  an  integral  equation  expressing  the  continuity  of 
the  tangential  magnetic  and  electric  fields  across  a typical  slot  in  a unit  cell. 

The  kernel  of  the  integral  equation  in  a sum  of  two  unit  cell  dyadic  Green's 
functions  G^,  and  ^ . The  former  (latter)  yields  the  magnetic  field  due  to 
an  elemental  magnetic  current  placed  in  the  slot  location  in  the  exterior 
(the  interior)  of  the  unperforated  structure, in  a unit  cell.  That  is  to  say, 

_Gg(r , i^P)  is  the  unit  cell  dyadic  exterior  G.F.  for  the  unperforated,  perfectly 
conducting  cylinder  and  Qj  is  such  a G.F.for  the  interior  unperforated  coaxial 
cable.  Both  G_(r,  /)  and  £T(r,  r')  satisfy  the  Floquet  conditions  which  implicitly 
include  all  interactions  between  neighboring  slots  : 


GE(r+xod,  _r#;  g)  = e ^dgE(r,  r*;  8) 

SI(£+£0d>  L<  8)  = e'jBd  £x  (*.  £,  3 ) • 


(1) 


The  integral  equation  reads 

JT^E 8>  + £l(r’  r’  0il  ‘ -X^o  dS>  = 0 ’ L in  the  slot  , (2) 

slot 
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where  E(r<)  is  the  unknown  slot  electric  field  and  p is  the 
radial  unit  vector  at  r ' . 


Application  of  Galerkin's  procedure,  i.  e.  an  expansion  of  the 
slot  field  in  terms  of  a linearly  independent,  preferably  orthnormal, 
finite  set  of  basis  functions,  and  the  requirement  that  the  resultant 
vector  function  (the  integral)  be  orthogonal  to  each  member  of  the  chosen 
basis,  yields  a homogeneous  linear  system  of  equations  for  the  unknown 
expansion  coefficients.  The  dispersion  relation  is  obtained  by  setting  the 
system  determinant  to  zero.  In  particular,  assuming  only  a sinusoidal 
distribution  e in  the  slot,  one  finds  an  approximate  dispersion  relation 


J J l0(±)x £0  • [£E(i’-/;  B)  +S?i(-’-/;  B >]'lo(-,)x-4  dS  dS'  = 0 

slot  slot 


(3) 


which  is  to  be  solved  for  the  unknown  propagation  constant  Pas  a function 
of  frequency.  It  is  seen  that  the  basic  ingredients  in  (3)  are  and  Qj 
The  expression  for  Gj  are  obtained  in  a standard  fashion  along  the  lines  of 
Appendix  A.  The  major  problem  in  this  study  is  the  determination  of  G,-,  . 


•TJOKSPSW'': 


3.  The  Exterior  Unit  Cell  Dyadic  Green's  Function  Gr(r,  r^P) 

— xl,  — — 

For  the  unperforated  circular  cylinder  the  unit  cell  exterior 
dyadic  Green's  function  is  represented  as  a sum  of  two  contri- 
butions 

£E  = = EF  + £EG  (4 

where  Grr  is  the  unit  cell  exterior  dyadic  Green's  Function  for  the 

— r_«  -T 

unperforated  cylinder  in  free  space  and  G^,^  represents  the  effect  of 

the  ground  and  the  multiple  interactions  between  the  cylinder  and  lossy  earth 

interface. 

An  approximation  for  G„_,  is  now  made  which  considers  only  a 

— F.  G 

single  ground  reflection.  This  approximation  should  be  adequate  , except 
when  the  cable  is  very  close  to  the  ground. 


3a)  The  G£F 

As  shown  in  Appendix  A,  is  derived  from  two  axial  Hertz 

potentials  TT<  and  TT“  which, for  the  slow  wave  case  to  be  discussed  in  this 
report, are  of  the  form 

a;  K (,  *■)"*  e',e"X 

n,  m=-  « 

n"(r,  r';  8)  = V A*  K (t  p).'*™  ” <x<d/2 

n,  m - -oo 

where  & - 04  n and  T = j 0^  - are  the  space  harmonic  axial  and 
n d n y n 

transverse  propagation  constants. 
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In  view  of  Floquet's  Theorem,  the  validity  of  (5)  extends  also 
for  -«  < x <«  . Thus,  a typical  term  in  (5)  represents  an  axial  electric 
or  magnetic  Hertz  potential 


A K (Tp)e-jmCPe-^X 
m m K ' 


(6) 


of  a multipole  slow-wave  (3>k)  phased  line  source  located  on  the  cylinder 
axis.  In  view  of  the  linearity  of  Maxwell's  equations  it  suffices  to  solve 
for  the  disturbance  due  to  lossy  ground  in  the  case  of  a monopole,  i.  e. 
for  an  isotropic  electric  and  magnetic  current  phased  line  source.  This 
is  because 


K (T  p)e 
m r 


-jmcp 


-ip?)]m[K0(Tp)]=LIt'Ko(Tp)  - 


Therefore , the  expressions  for  the  fields  excited  by  an  mth  multipole 

line  source  above  ground  may  be  generated  by  the  application  of  the  linear 

differential  operator  (L)m  to  the  fields  due  to  an  appropriate 

isotropic  phased  line  source  above,  and  parallel  to, a lossy  ground  interface. 

3b)  Basic  Integrals  for  Q . 

In  Appendix  B it  is  shown  that  the  fields  of  an  electric  or  a magnetic 
line  source  above,  and  parallel  to,  a dielectric  half  space  may  be  derived 
from  two  Hertz  potentials.  From  (B  14,15,26  and  27)  one  observes  that  the 
basic  types  of  integrals  to  be  evaluated  in  order  to  determine  Gr  are 


*2  = J k e -e  -t)  , lm^2<0 


l* 

*2  = J 


- jTJy  -JH^z+z') 
e e 

(K2+G  Hl)  *1 


dn 


(«) 


Other  integrals  that  appear  are  obtained  via  the  operations  or 

applied  to  1^  or  I^.  As  indicated,  we  shall  restrict  our  considerations 

to  the  slow  wave  case.i.e.  , such  that  for  e = kyec  (ec  being  the  dielectric 

constant  loading  of  the  coaxial  cable),  the  spacing  d/X  is  sufficiently  small 

so  that  for  all  n,  ^ >k2  , At  higher  frequencies  the  cable  will  become 
n 

leaky  when  one  of  the  space  harmonics  becomes  fast.  The  expressions  for 

a leaky  cable  may  be  derived  via  minor  modifications  of  the  results  given  in 

this  report.  I.  and  1?  maybe  further  simplified  as  follows: 

1 ^ 


03 


kZ(e-l) 


r 

J 


V9 

*i 


-iny 

e 


-j*..  (z+z') 
e dn 


jk2(e-l) 


itz-  ^ */y2+(z+z')2) 


kZ(e-l).Je 


2 -jny 

— e e 


jHj  Z 


dn 


where  Z = z + z and  is  the  Hankel  function  of  the  zeroth  order  and 

o 

2nd  kind. 

The  integral  may  be  rewritten  as: 


Z2  = 


1 

7T 


(»,-€«,)  -i'15'  ->*iz 

' 2 1 e e , 


~2 — r 

n -a 


dn 


(9 


(f 
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where  a - k ^ - - p 


^ k jj—  - B ; Im  a < 0 


The  poles  at  Ti  = ±a  are  the  so-called  Sommerfeld  poles. 
P 

The  basic  integrals  to  be  evaluated  are  therefore 


F<v.  Z)  = J *2  Z-z 


-j'ny-j^z 


“ -jTiy-i^z 

G(y,  Z)  = j S 2 dn 

-»  ft  'a  )*i 


r *2  e 


-jTly-j^Z 


H<Viz>  = J -z-r dTi  • 

_»  (T-a  ) Kj 

F(y,  z)  will  be  evaluated  as  an  asymptotic  series  in  terms  of  the  canonical 


integral 


jTly-jKjZ 


dri  = n k2-P2  J y2+z' 


This  expansion  is  valid  for  sufficiently  large  distances  and  will  be  de- 
scribed in  Section  4 . 

The  cano  lical  integral  G(y,  Z)  is  evaluated  exactly,  and  H(y,Z)  is 
in  turn  expanded  in  an  asymptotic  series  in  terms  of  G(y.Z)  and  its  de- 
rivatives. The  integral  F(y,  Z)  represents  the  TE  mode  contribution  , 
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while  G and  H , which  contain  the  Sommerfeld  poles,  represent  the 
TM  mode  contributions  (with  respect  to  z). 

4)  Image  Series  Representation  for  the  TE  Contribution  F(y,  Z) 

We  are  considering  here  the  slow  wave  case,8>k  and  define 
now  the  following  quantities: 


r~2 — z rz 

\ Vp  -k  51  °-  *i  = -WY1 


-jV  , 


Re  y>o  , 


(17) 


7 kVl>  , >2  -VV22  ^ 


, I m K ^ < 0 


(19) 


k J €-1  , I m < 0 


(19) 


K(y,  Z)  may  now  be  written: 


P e'iT|ye'Y7 
F(y,Z)  - jj  k2  £ -£ dT| 


(20) 


When  k7 -1  the  integral  in  (20)  represents  the  modified  Bessel  function 
,r_2 2 - Y Z 

2 K (vjP),  u - y y +Z  . In  view  of  the  factor  e , the  major  contribution 
to  the  integral  arises  from  the  neighborhood  of  Tj-0.  As  a result,  when  | e| 
is  large,  ^ Y^-"'2  in  (20)  may  be  expanded  in  a Taylor  series  about  Y “0: 

, . . j , I m Y'2  < 0 


71  71 


Vy2-V  Y 2 L*  ~ 17 


1 Y2  Y4 


Hyz 


(21) 


and  the  resultant  series  is  integrated  term  by  term.  The  integration  is 

2 

facilitated  by  the  presence  of  powers  of  Y . whereby  Y — y acting  on  the 

bZ^ 

basic  integral. 
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M'-fl  -i  »%, ** 


Thus , 

F(y,z)  2Y2  LKo(YiD)  - ^2- 


az 


(22) 


J? 


ft 


Jtt* 


The  series  in  (22)  is  a multipole  expansion,  the  first  term  representing 
the  monopole,  the  next  a quadrupole  etc. 

As  a result, 


2 2 
Y,  Z 


, r Y1  Z , Y1  2Y1Z  V 

F(y , Z)  ~ 2Y2  {k  (YjD)-  -yLW*  2 " K1(YM~ 5 J. J ' 

2Y~  c P 


2Y' 


2„2  „4^4 


2„4 


8Y 


, r 3Yi  24Yiz  Yiz  24YiZ  ^ 

HKo<YA“2 4 + — ^-  + —5— i 


r 6Y.  6y,3Z2  48YiZ"  8y^J 77  48y1Z‘1, 

+■  ki(vic)L — T 3 5 + 5 + 7 jj 


3„4 


+ . . . 


/ VlW,  y2  (Yl^Kl(flp)  A 2v2  \ 

d(yip)-W  a-tt’  tt— v--z-y 


”2Y2  \K 


2 '2 


Z(Y2°) 


2 2 2 


K 


o'YK-t:)  [<V>‘2  (t-’Vv1- 7-) 4 ttVtj-)  \l-y  J 


-2, 


-2  , Y.  2 2 

3 , 1 , „ y 


■ 2 2 


-(YjoJKjfYjC)  (Y2p)  lT'(Y2p)  6(Y2r) 


-q 


2 2 2 - 2 ^ ^ Yp  - 

+ (-^-)  (i--V)  + 6<v)  (1"V)  ]} + °l(_y7)  (Y}  j 

' 9 r 2 


p+q  - 6 
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(23) 


(24) 


\ •« 


A 


» 


When  the  cable  is  not  too  close  to  the  ground  surface,  then,  for  observation 

points  on  the  cable  (or  on  a typical  slot)  y/,  <<1,  and  one  may  neglect  in 

?. 

(24)  all  powers  of  ( -f-)  and  higher,  or  retain  only  the  y independent 
terms  plus  Of  — I . 

V ) 


For  y=0  one  has 


F(0,Z)~2Y2lKo(Y1Z)-Ko(Y1Z)T-(-)- 


1 , V (YjZ)Kj(YjZ) 


2 2 (Y  2 Z) 


Y,  2n  Y,  p q 


K (Y,Z)  Y,  \ o 1 Y1  N.  1P 

A-V-(t-)  +(V1Z)K1(Y1Z)(Y2Z)  )V-(V2Z)  ^(—^(-H^z 


:2V2{k(ViZ){. 


Ko(YlZ)  Y1  (YjZJK^YjZ)-! 

5 ( "v  ) + 5 i 


-2  2 

1 Y,  Y,  p -q 


2 ' Y 


2 2(Y2Z) 


JH/  t 1 '1  1 r ^ 

—J\1+_r(Y2z)  + t^;)+0<t->  (y2z)  p+q  6 • 


It  is  seen  that  the  asymptotic  expansion  (25)  contains  two  small  parameters 

Y1  2 -2 
( y — ) and  (Y2  Z) 


Example: 

For  je-1 1 = 10,  = 1.6,  Z = z+z'  = 0.  2\ 


Y,  Z = l.  57,  | Y,Z|  = 0.  06 


1 1-1 

~Z~  Y -> 


= 0.  07  . 


He  nc  e , 


K (Y,Z)  Y,  (Y,Z)K  (Y  Z) 

— T - ( Yi)  + — -T— 

Y2  2(Y2Z) 

K^^Zl 


If  one  neglects  the  3rd  order  term. when  its  magnitude  is  say  10% 
of  the  second  order  term,  i.  e.  , when 

-2  . Y.  Z 2 

4|v2z>  tr(Tpr)  'i  °-'  ■ 

then  for  |c-l|  =10  and  = 1.  6 , 

/ , 

.-+  z . must  exceed  0.16.  This  is  the  case  when,  for  z=z  , 

Jt  >0.08,  i.  e.  , the  height  of  the  cable  should  be  greater  than  approx.  0.  1' 

Ihus, for  given  0/k  — of  the  unperforated  coaxial  line  and  for  given 

VP 

electrical  properties  of  the  lossy  ground, the  validity  of  the  asymptotic 
expansion  of  F(y,  Z)  (24)  or  (25)  is  limited  to  distances  z not  too  close 
to  the  ground.  Formulae  (24)  or  (25)  permit, in  each  case, to  assess  the 
validity  of  the  expansion  i.e.  , of  the  minimum  distance  of  the  line  source 
(cable  axis  to  ground)  for  a given  value  of  and  6 . 

5)  Evaluation  of  the  TM  Contribution  G(y,  Z) 

Since, for  the  TM  polarization,  the  air-lossy  ground  interface 
supports  a Sommerfeld  surface  wave,  it  is  clear  that  an  image  repre- 
sentation of  G(y,  Z)  will  not  be  useful. 

The refor e, one  has  to  resort  to  a type  of  evaluation  which 
explicitly  exhibits  the  surface  wave  contribution.  Again  we  consider 
here  only  the  slow  wave  case,  for  which 


yY1  +y  Z - )T|y 

- jG(y,  Z)  j -2  2 ---  dr\  = j % 2"  7F A (1~ 1117 ' 

-»  (T)'  -a  jy^+ri2  (11  -a  )./^4Tf 


- /"Z 2” 

VY1  +T1  Z 


2 2 <26> 
-V-+...)dn. 


Thus,  since  the  odd  powers  of  y do  not  contribute 


G(0,  Z)]y2+0(y4}  , 


(27) 


r Yl. 


G (y , Z)  = G (0,  Z)  +L-fG(0,  Z)- 


SZ 


where 


G(0,  Z)  G(Z) 


■fil 


+T)  Z 


7.7  2 

+ T) 


dri  . 


(28) 


It  is  therefore  sufficient  to  evaluate  G(Z),  whereupon  G(y,  Z)  follows  via 

(27). 

For  the  purpose  of  contour  deformation  that  follows,  it  is  necessary  to 

discuss  the  properties  of  singularities  in  the  complex  T|  plane  of  the  integrand  in  (28). 


5a ) The  Top  Sheet  of  ^ > 

2 7 

in  the  Complex  ri  Plane. 

, W e have 

2 , 2 e 
a =k  yy-jt-  ■ 

1+t 

j 

- p , lma<0  (tor  definiteness)  . 

(29) 

: Re  J v i +fi2 

! 

> 0 , on  the  entire  top  sheet 

(30) 

g = 3r+  jBi 

> 8 > 0,  £L  < 0 ( for  propagation  in  the  +x  direction 

(31) 

with  a simultaneous  attenuation  in  the 

direction  of  propagation  ) 

e = e +je.  , 

r i 

e >0  , e.  < 0 
r l 

(32) 

rjr,  2 [2 — 7 :2'  : 

l ~ 'J  “ ” k = V"r  "Bi  *k  +2]°rei  • (33) 


Si  nee,  on  the  top  sheet  Re  Y ^ - (Re  J > 0 , 
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A 


•.TiiS#s2P  WSK&2&  * 


f 1 

> 

I < 

: V. 


I 2 T . 

Along  the  branch  cuts  Y^+T|  is  purely  imaginary.  To  determine 
the  proper  sign  on  either  side  of  the  lower  branch  cut  we  observe 


i~2 — 2 

that,  on  the  top  sheet,  the  sign  of  Im^yv^  +T)  on  the  right  bank  of  the 


lower  branch-cut  is  that  of  Im  at  T|  - 0 i.  e.  of  ImY,  which  is 

negative  as  shown  previously.  The  sign  of  I m +T)  changes  across  the 

l~Z  T 

cut  and  therefore  ^Y^  +T1  becomes  positive  imaginary  on  the  left  side  of 
the  lower  cut,  as  indicated  on  Fig.  4 . 

With  the  above  c onside rations, we  deform  the  contour  in  (28)  from  the 
real  axis  around  the  branch  cut  and  capture  the  pole  at  T)  ^a.  The  de- 


T~2 


’ 2 2 

formation  is  permissible  because  on  the  entire  top  sheet  Re^j  y^+r\  >0  . 

(Actually,  since  the  term  e ,T1'  has  been  removed  from  the  integrand 
in  G,  the  contour  could  be  just  as  well  deformed  about  the  upper  branch 
cut  in  Fig.  4. 

The  result  of  the  contour  deformation  is: 


G(Z)  = - 2tt  j 


-R 

. e 


■R 


+ T)2  Z 


J 


n 7 " 2 2 r~2  2 

-a  )Vvi+T» 


dn 


As  shown  in  Appendix  C,  the  poles  at  Tj  = ±a  are  located  on  the  proper 

j~2 — 2 

sheet,  i.  e.  Re^Yj  ta  > 0 . 

Changing  variable  via 

TJ'J.  r . CdC  CdC 

VV  ^ = C , i.  e.  dn  = = 


"2  2 
V C -Y 


nz~r 


and  remembering  that  ImT)<0  on  the  cut,  so  thatlm^1*  - Yj  <0  on 


the 


contour  C,  , one  has 

b 


(34) 


(35) 


T*  * 

* 

* 
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m * .»•# 


v/'» 


- 1® 


G(Z) 


-JV"  Z 

jTT  e 


e"jtZ  dC 


~T~ 2 ZT 


Yj  +a 


U -Yra  ),/c-.y 


7T~2 


I m C - Yj  <0 


Upon  substitution  C = j i 


r~z — T 

(°rVvi+T1  i*) 


one  finds 


z 


G(Z) 


-jne 


•Ay  + a 


* 3 


-jCZ  ,, 

e d>o 


2 2 Z , / .-2 — -Tf 
(--  rYx+a  )^-(/+y2) 


-JTT  e 
a 


A 'i 


r.  +a  Z 


e'^Z  dS 


"Z — F 


.-2  2 


wa 


) V 


7Z- F 


“2 — Z 

Re  ^ Y,+a  > 0 


~2 — T ~2 — Z 

Since  ; = -yY.  +T)  , the  proper  sheet  Re,yv  ^T)  >0  maps  into  the 


lower  half  5 plane  I ml  - 0.  The  upper  half  ? plane  corresponds  to  the 


T~T 


improper  sheet  of  the  ^ ^ +tj  . Let  the  pole  of  the  integrand  at  : ^ 

/"Z"”2 


have  Im£  <0  i.  e.  let  Re^'Yj  +a  > 0 . The  other  pole  is  then 

p2 


, ~2 F 

j^Yj+a  with 


,r'2— Z 


ReV"-i +a  <0  ■ 


In  this  fashion  the  resulting  contour  of  integration  has  been  again 
brought  to  the  real  axis,  and  therefore  the  integral,  which  is  of  ar 
inverse  Fourier  transform  variety,  will  be  further  simplified  lay  use  of 
convolution 
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mrJF9ZrX*F'2V 


CO 


j u(Z-0  v(C)  dC 


1 r -i?Z  _ 

if  J U(S)  V(3)e  d5 


u(C)  - ^ j U(S)e"jib  d§ 


v(C)  = j V(§)e'jCS  d§ 


Therefore,  from  the  integral  in  (38)  one  has 


u(Z-fi)  = 


i p e-jS(z-Q 


r+Yj  + a‘ 


„ i p e-i^ 

V(C)  = Tu  J “7= 


To  evaluate  (41)  for  Z>G>  one  closes  the  contour  at  infinity  in  the  lower  half 

plane  and  takes  out  the  residue  at  5 , • For  Z<C  the  contour  is  deformed 

pi 

in  the  upper  half  plane  and  the  pole  5 ^ 's  captured  with  the  overall  uniformly 


valid  result 


u(Z-C)  = 


yf+a2  I z-cl 


Re  /Y,  +a  > 0 . 


Evaluation  of  v(£)  yields 

V(C)  = 4 Ko(Yi  !§l) ' (44 

where  K is  the  modified  Bessel  function  of  the  second  kind  and  zeroth  order, 
o 


A* 


Consequently  one  has, 


G(Z) 


- jTT  e 


r-y T _ f~ 2 T 

Z * 7Va 


V Yi 


I z-C I 


rz  2~ 

VYi+a 


K 


rz — r 

VYi+a 


0(Y1l;l)dc, 


R 


Pe./Vr  + a >0 


(4 


Re  Vj  >0 


Im  V.  - 0 


To  further  simplify  the  integral  I(Z)  in  the  (44)  , one  has 


r~Z — Z 

• ^/Y,  +a  Z Z 


I(Z) 


rz — z 

VVa 


r . 

j dC  e Ko(v|ri)d;+ 


^Yj+a  Z » '.,2,_2 


r -Vl+a  C 

1 c 


— 7 7 

VY+>  z 


«„('<!  It  lw. 

(4 


or 


Z — z 


+ a Z 


Z Y^  + a2  * 


I(Z)  = 


■7yT^ 


r_Z — Z - 

r“  _VYi+a  ^ V T 

LJ  e Ko(Y1C)d;+Je 

O r> 


K (Y.C)dt 

O X 


+aZ  Z 


rz — T 

VVa 


ot>  r~2 2 7 FZ T 

r , V Va  C ? ,/Va  S 

Lj  ‘ WC  -J  « Ko<Y>d;j  = 


2 cosh 


ff*? z * w c 


y 


2 T 

Yj  +a 


Ko(Y1C)dC  + 


>l+a  z z 


rz — z 

vva 


j e 

o 


-R 


^(Yi^dC- 


+a“  z z c 


J 


7~Z — Z 

VVa  ° 


Ko(Y1QdC  (4 
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But  from  l^] 


00  / 2 , 2 - 
(•  _ _y  Y 1 ^ a 10 
i d»e 
«J 


Ko(Yi;) 


WYa 


ReVVa  >0>-  ReY,  . 


where  trie  principal  value  of  fn  is  understood,  ReVy^+a  >0  , and  the 

r~7  Z 

convergence  of  the  integral  is  assured  since  both  Re./Yj’+a  >0  and 
Re  v > 0 . 

Thus,  the  evaluation  of  G(Z)  is  now  reduced  to  the  evaluation  of  the  tw  o 
finite  integrals 


WT+a  C 

d£  e~  KQ(Y1C)d;  . 


For  small  values  of  V^Z  i.  e.  , | >^Z  | < < | , Kq(w)  may  be  approximated  by 
4nV  w/2  where  Y is  the  Euler  constant.  This  leads  to  expressions  in- 
volving exponential  integrals.  Otherwise,  in  the  dispersion  relation, which 
in  view  of  (26),  (27)  and  (Cl),  contains  I and  their  derivatives  with  respect 

to  Z , the  most  rapidly  convergent  procedure  valid  for  any  value  Y^Z  is  to 

i 2 ~Z 

evaluate  numerically  with  Y^  = ,yp-k  as  a 3 dependent  variable,  and 

T Z - ik 

'■v  co  — > — as  a B independent  parameter. 


Vl  + e 


Altogether , 


G(Z)  = j ^ 


— j TT  e 


(Y^  +a  Z 


, , ,,/Y,  +a  Z, 

2cosh(V  1 ) 


\ T 


Vva< 


z f ^ 

-=ifT  Joe 

VYi 


K0(Y^)dC  + 


e^+a  z r -vV+a“  C 

i e K 

° 


With  G(Z)  given  by  (49),  H(Z)  is  calculated  in  a similar  fashion  as  F(Z), 
by  expanding  via  (21). 

Conclusions 

A rigorous  formalism  has  been  presented,  that  permits  setting  up 
of  the  dispersion  relation  for  a periodically  slotted  coaxial  cable  above, 
and  parallel  to,  a lossy  ground  in  the  surface  wave  region.  A simple 
modification  of  the  procedure  will  allow  to  account  for  one  or  more  fast 
space  harmonics  and  thustoobtain  expressions  valid  for  periodic  leaky 
cables.  The  expressions  may  be  easily  analytically  continued , (in  contrast 
to  those  in  Ref.  [2],  in  the  case  when  the  Sommerfeld  poles  as  a function 
of  complex  propagation  constant  3 would  be  required  to  cross  the  real  Raxis. 

An  interesting  observation  is  appropriate  here.  Consider  expression  (28  ), 
in  the  complex  R plane  with  its  Sommerfeld  poles  at  R^±a  in  the  Im  a < o 
which  are  a function  of  |3  . As  recognized  in  Ref.  L2l  , the  poles  move  in 
the  complex  T)  plane  as  a function  of  the  complex  axial  propagation  constant 
8,  which  is  governed  by  the  dispersion  relation  of  the  periodically  slotted 
structure.  If  the  pole  at  Rp  = a,  as  a function  of  3,  crosses  the  real  axis  into 
the  upper  half  plane,  the  pole  at  R = -a  crosses  the  real  Raxis  simultaneously 
into  the  lower  half  plane  and  the  representation  (28)  is  a discontinuous  function 
of  3(k).  Such  functions  however  are  not  admissible  in  the  dispersion  relation 
and  therefore  the  function  G when  R^=  a(g)  approaches  the  real  axis  from  below, 
must  be  analytically  continued,  as  a function  of  g . A similar  situation  occurs 
in  Landau  damping  L5],  and  the  procedure  is  analogous  here.  As  the  poles 


± a cross  the  R.  = 0 axis  the  contour  is  indented  as  shxvn  in  Fig.  5a.  After  the  pole 
R = a crossed  into  the  upper  half  plane  the  appropriate  contour  is  shown  in 
Fig.  5b. 


(a)  CONTOUR  FOR  SOMMERFELD  (b)cONTOUR  FOR  SOMMERFELD 
POLES  JUST  CROSSING  POLES  AFTER  CROSSING 

THE  REAL  AXIS.  THE  REAL  AXIS  (FROM  BELOW). 


Fig.  5.  Analytic  Continuation  of  G(y,Z,3)  . 


The  deformation  of  contour  provides  a simple  means  of  analytic  con- 
tinuation of  G(y,  Z,  3)  whereby  the  appropriate  pole  'H^-a  is  tracked  even 
if  it  crosses  into  the  upper  half  T plane.  Physically  .after  the  pole  f]  = a 
crossed  the  real  t|  axis, so  that  now  I mp ^ > 0,  such  a pole  represents  a leaky 
wave  in  the  y direction.  This  is  obvious, when  in  (26)  the  contour  would  be 
deformed  into  the  lower  branch  cut,  prior  to  Taylor  series  expansion  of 


the  e 


References 


1.  J.  R.  Wait:  Electromagnetic  Field  Analysis  for  a Coaxial  Cable 

into  Periodic  Slots. 

2.  R.G.  Olsen  and  D.  C.  Chang:  Electromagnetic  Characteristics  of  a 
Horizontal  Wire  Above  a Dissipative  Earth,  Parti.  Scientific 
Report  No.  3 Electromagnetics  Laboratory  Dept,  of  EE,  University 
of  Colorado,  October  1973. 

3.  M.  H.  Chen:  Electromagnetic  Wave  Propagation  Characteristics  of 

Periodic  Structures  with  Screw  Symmetry*  Ph.  D.  Dissertation, 
Polytechnic  Institute  of  Brooklyn,  June  1969. 

4.  W.  Magnus,  F.  Oberhettinger,  R.  P.  Soni:  Formulas  and  Theorems 

for  the  Special  Functions  of  Mathematical  Physics,  p.  449,  Springer  1966. 

5.  B.W.  Roos:  Analytic  Functions  and  Distributions  in  Physics  and 

Engineering,  Wiley  1969. 


?J«0r  3Z*F'Zjr  ..  v) 


Appendix  A 

Summary  of  Radial  Line  Formalism  for  a Unit  Cell  in  Cylindrical  Geometry 
(based  on  Ref,  13  ~!).  

p When  a scalar  and  a vector  multiplic  ation  with  a radial  unit  vector 

is  applied  to  the  Maxwell's  equations 
VxH  ~ E 

VxE  = -jxp  H-M,  J 

— o—  — 

the  following  set  of  coupled  equations  arises  I'm-  the  transverse  to  t 
field  components  E^  and  . 


v 7 

-vpEt=  • -Mt  x -O 


?(Ht  j'lje  . 1_  + 


(.  x E.)  , 
— o — t 


, „ X 10  6 

w here  7 x — — — r z --r — 

t o p 0 y O OZ 


The  unit  cell  orthonormal  vector  mode  functions  that  scalarize  this  set  oi 
equations  are 


e'  (r)  = - -L—  » + z 

£mn-  y^rr-d  L V1  -S 


pm  -jm$  -j8  z 

n • n 

§ 4-  z e e 


h'  (r)  — I * ■ e e 

-mn  - JZWT  w D 


■ jm$  -jj3  z 


E-type  modes 

(H  -0) 

z 


, r , - im$  - jS  z 

e"  (r)  = j — * ! e e n 

-mn"  c/ZttcT~  L P J 


h(r) 


jmi  - j 6 z I z 


H-type  modes 

(E  =0) 


i r r m - -pris  - 1 s 

_L_  | - _J?_  « + z . e e n 
rnr  - *2e 


• x**<*»~  -*  **  »<• 


* '+i 


From  (A  6)  one  finds  that  for  M =0  v' (p)  and  I ^n(p)  satisfy 

the  Bessel  differential  equation 


mn 


fid-  d 

Lt  3TT  (p  SZ 


2 m 

— i / t r\ 

dp  n 


V (P)> 
mn' 


■mnl^ 


= 0 , 


A(10) 


while  the  generator  voltages  v'  and  v"  are  obtained  for  the  case  of 
& mn  mn 


• o / ' 7 7/  " \ 

transverse  magnetic  current  sources  via  an  expansion  Mx£q  = h 'Vmn-mn+Viiin-mii' 

so  that  using  (A4) 


v*  (p)  = I i M ( r )•  h'  (p)  ds 

mn'  ’ J J — t — — mn'K/ 

unit  cell 

cross-section 

A tangential  magnetic  point  source  M^tr ')  on  an  exterior  of  an  unperforated 
coaxial  surface  in  a unit  cell  may  be  represented  by 


A(ll) 


M,(r  ) = u M 
— t — o 


Mr  -a)  tj'C-V  ) f)(z- /■  ) 


A(12) 


where  u is  a unit  vector  in  the  direction  of  M,  , 
o — t 


Hence, 

Vmn(fl)  = M %*  — m rS—  ^ 6(o'a)  = Vmn  ^'a) 
v"mnW  " Muo  ‘ — mn^— £<^a>  = Vmn  ^P'a) 


A(13) 


Upon  integration  of  equations(Ab)  one  finds  that 


_(V#  (a)-V'  (a))  v' 

mn  mn  mn 


-(V  (a  )-V  (a))  = v 

mn  mn  mn 


A(14) 
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and  since  V*  (a)  V'  (a)  0,  in  view  of  the  perfectly  conducting  cylindrical 

mn  m n 

surface  at  p a, 


V (a)  -v 
rnn  mn 


V’  (a)  = - v 
mn  mn 


A ( 1 5 ) 


In  vieu.  of  (A10),  for  the  exterior  of  a cable  in  free  space 


V'  (p)  = A H(2)  (n  p' 
mn  m n 


A(16) 


Hence, 


V'  (a+) 
mn 


A H (Z* (h  a) 
m n 


A(17) 


and  the  refore, 


V (p) 
mn 


Mu  • h*  (r  ) 
o — m n — 


H ' ’ (k  a) 
m n 


h(2)(^  p) 
m n 


A (1H ) 


Thereto  re  .from  ( ) 


-Mu  • h ' (r  ) H ' (k  c) 
o — mn  — m ' n 


I (p)  - 


H ' ' (a  a) 
m n 


H(2)(ht)  A(19) 

. . ...  m n 

„ .hw  (r)  . 

n m n 


To  calculate  l"  ( :)  one  has  from  (A10) 
mn 


I"  (r)  - R H(2)  (h  r)  , 
mn  m n 


A(20) 


so  that 


i"  (a+)  = B H (2)  (h  a) 
mn  m n 


A(Z1) 


m*  v. 


and  therefore 


H(2)  (k  p) 
i"  (p)  I"  (a  + ) - 

mn  H{2)(*a) 

o n 


A(22) 


From  (A6)  one  has  via  (A  15) 


//  1 d T"  , N 

V (r\  i T—  I (p) 

™(  ) (r)  y"  (p)  dp  mn 

mn' mn 


h<2>(*  p) 

(a4)  ™ 

A mn'  a\ 

n H m ^nd 


H(2)(n  a) 

V"  (a+pJ^-l"  (a  ) ^ 

mn'  Kn  mn  h(2>(k  a) 

m n 


-Mu  • h (r  ) 
o — — 


A ( 2 3 ) 


Hence>  „ / (2), 

//  , +1  MV  — m n ^ Hm  Va) 

rmn(a  ) = ‘ j iX1  pa  a.  "n  H(2)('K 

m n 


A(24) 


I (p) 

mn 


M • h“  (r  ) 
— t — mn 

jx  p-L  a 


h(2)(V) 

m n 
m n 


A(25) 


Consequently, 


H.  (r , r ) = G(r,  /)*M  - 
-t' -FF 


, Z H^2^  (a  p)  .* 

I S Y ' ("'l  h'  (r)  h (r  ) 

U n p h(2)  (h  I7“mn  ~mn~ 

mn  = -ai  m n 


H<2)(h  0) 


1 


L-J  11  /I.  v 

mn=  -oo  (Kna) 


A(26) 


.#*•■(*  .i  «*..(»*  *■ 


Hence, the  transverse  part  of  the  exterior  dyadic  Green's  Function  for 
a coaxial  cable  in  free  space  is  given  by 


H(2>(HP) 

Jr.  — — mn^r^  + 


m n - c 


Hv  '(k  a) 
m n 


" H(2)(k  r) 

\ Y"  (a)  ™ " h"  (r)h'T(rV 

n H 2 ( k a)  ~mn  - 

in=-co  m n 


vv  ith 


£ P 

n k 


a 

Y"  n 


n jtu  pp 


It  can  be  readily  verified  that  G can  be  derived  from  two  scalar 


=EF 


potentials  TT  and  tt 


n 


' = i,jje  y 


, H^(r  p)  -jm$  -jg  z 
1 m n n 


«/2nd  l"  ^n^  H ^ 2 \ x a) 

mn--®  m n 


, , . , H^2\n  p)  - jm § -jg  z 

1 1 \ 1 m nK  ’ n 

TTTT / T. TTV e e 


* 

i? 

t 

I 

V 


TT  = - 


j^P  yzffd-  ^ Kn  1^[k  a) 
mn=  - oo  m n 


Thus,  both  tt  and  tt  consist  of  terms  of  the  form 

r „(2),  . ’jenZ 

C H (k  p)  e e , 

mn  m n 


§ 


each  representing  a Hertz  potential  of  an  electric  (TT  ' ) or  a magnetic 

m n 


// 

(nmn)  phased  line  source.  Expression  for  the  interior  dyadic  Green's 


Function  Gj(r,r  ;g)  may  be  derived  in  a similar  manner. 
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A(29) 
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Appendix  B 


Hertz  Potentials  for  Phased  Line  Sources  Above  a Dielectric  Half  Space. 


The  potentials  for  an  electric  line  source  appear  in  [2]  and  were  originally 
obtained  by  J.  R.  Wait.  For  completeness  we  rederive  them  for  the  electric 
line  source,  and,  in  addition  obtain  expression  for  a magnetic  line  source. 

1.  Electric  Current  Source  at  _P ' = (y\z')  above  the  interface. 

We  assume  that  E and  H are  derived  from  two  Hertz  potentials  II 'xo  and 
Il"x  such  that 


— o 

In  the  air  region 


In  the  dielectric 


H ,i_  = i uce  VxII /+x 
— J o — o 

E ' + = VxVxII , +x 


— o 


_//-+■  . _ n //  + 

E = - ] ULU  x 

— J o — o 


h"+  = vxvxn"+x 


v2n,+  + k2n,+ 


iuue 


v2n"+  + k2n//+  = o 


H 7 ’ = iuuee  Vxll 7 x 
J o o 

E ' ' = VxVx n 7 x 

— — o 

E//_  = - j uup.  VxlT  x 

H"~  = VxVxlT'x 

— — o 

v2n'_  + k2en'"  = o 
v2n/,_  + k2en"‘  = o 


B ( 1 ) 


The  primary  electric  line  source  fields  are  generated  by  11^  . The 

scattered  fields  in  the  air  space  are  derived  from  Il'+  and  Il"  + such  that 


2„/ 


n,+  = n/  + n'+,  n"+  = n",  with  v2n/  + k^n: 

inc  s’  s me  inc  juue 

as  well  as  V2n  "*"  + k II  =0. 

s s 

In  the  dielectric  half  space  one  has  II  = IIe 


^ and  V2n'+  + k2n'+  = 0 
s s 


o 


both  of  which 


satisfy  the  homogenous  Helmholtz  equation. 
If  2 = Joe"^0X6(y-y  ')6(z-z  ')xq 


B(2) 


then 
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A'tt.r 'X*r 


l -J9x 
oe 

j“eo 

4 h(Z)  ( 

4 o 

!'  2 R2  . /., 
yk  -9  V(y-y  ) 

1 / i 

-I  -j3x 

00  -j^y 

?«e  e 

-JHj |Z-Z  | 

oe 

UUf!  4lt 

J 

H1 

o 

00 

,/,2  e2  2 

= vk  - s - n 

lmn<0 


Setting,  in  accordance  with  the  radiation  condition 

, -i3v  r - j"ny  +jH,(z-z') 

II.  = e J A{T))e  e dr|  ; 0 < z < z 

inc  J ' 


ns+  = e'j3X  J (BW-  Mr\) 


-jv  y -jfi^z+z') 


e e 


n»+  -j3x  r , •J,|y 

n = e J C(n)e  e 

3 W 


-jTiy  -jK1(z+z/) 


n"  = e 

s 


,Sx  r -j^y  -j*?2 

= e J x J D(ri)e  e dr) 


-i3x  T 'jTly  'j*2z 

n"  = e J3X  J E(V|)e  e d dr) 


and  imposing  continuity  of  E^,  E^,  and  at  z = 0 one  obtains,  in 
view  of  (Bl)  and  (B4)  the  following  set  of  linear  inhomogenous  equations 
for  the  determination  of  the  transforms  B,  C,  D and  E: 

22  ? 2 “3*1 z 

(ke  - S^D  = (k^  - S^Be  1 B| 


22  ? 2 z 
(k  e - 9 )E  = (k  - 8 )Ce 


• Tjasra?* • 


-jH.Z  -j^jZ 

■3r|Be  -uu^HjCe  = -3r|D  + uu|Joh2E 


-J*,Z  -JH.Z  -JH.Z 

■ 2iii£  K,e  A + uue  >t.  e B- Brie  C = -U)ee  h0D  - BtiE  . B(8) 

o 1 o 1 o l 


From  (B3)  one  has 


A™  ° 4 ~t  », 
o 1 


Elimination  of  D and  E yields 


2 2 2 
2uJeo(e  - (l>  )A  = -^otC  • ' vl)  K + CKz{1  - (I)  JB  + BTl(e  ' 1)C 


B(  1 0) 


0 = 8ri(e  - )B  + uuia^  Ve  " \k/  /h1+h2\'V17/  )j 


After  certain  amount  of  manipulations,  the  determinant  of  the  system  A 
may  be  shown  to  reduce  to 


A - - V e - l — 


1 -(4)  AK1  + H2)  WChi) 


B ( 1 1 ) 


Consequently,  from  (B9),  (BIO)  and  ( B 1 1 ) one  has 


E - ° r 1 W j 1 

4n“£0  L Vk2  K2+exl  J , f3\2 

1 ■ M</ 


B(  1 2) 


I 3/'k  r 
o 1 


2TTk^-©  ) 


r i i i ji 

■ Lk2  + h1  H2+eH1  J H1 


B(13) 


+ = ' °-  e‘j0X  iH(2)/k2-02V(y-y,)2+(z-z')2  - H(2)(42-BW(y-y,)2+(Z+z') 
4&e  o 


) - 


f r i 11  -jTly  -j*i(z+z,)  . i 
J Lkj  + K^  " K.2+eKl  J e e 11  ) 


B(: 


The  expressions  differ  by  a factor  of  2 in  front  of  the  | j bracket  and  by 


factor  of  — r— in  front  of  the  integral  from  that  in  Eq.  (1.  2)  of  Ref.  [2]. 
-J  i 


T//  + o 


I S/k  e'J 


- i0x  ® 


) - 


r r_j_ . i in 

27  J H2+eKi  JXj 


-jly  jK^zfz') 


B( 


This  expression  checks  with  (1.4)  of  Ref.  [2], 
2.  Magnetic  Current  Line  Source. 


J = M e”^0x  6(y-y')  b(z-z')  x 
— m — o 


B( 


E and  H due  to  (Bl6)  in  free  space  is  derived  from  _xQ 
where 


v2n"  + k2n"  = , 

inc  inc  jiDp 


B( 


with  the  radiation  condition  at  infinity. 
Thus, 
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+*  -v. 


‘■'“jatrs R*T:'-Zf  k 


S U-  Wk -r  Vy-y')' +(z-z') 

JUUH  \4  o \ 7 7 


„ -jsx  : -j»i u-^' t 

Me  J e dr) 

4niuu  J k, 

o i 

«.  00 


A magnetic  phased  line  source  parallel  to  an  air  dielectric  interface 
excites  fields  that  are  derivable,  similarly  to  the  case  of  an  electric  line 
source,  from  two  Hertz  potentials  via  the  first  four  relations  in  B(l),  the 
difference  being  that  now  in  addition  to  (Bl7) 


n"+  = n?  + n"+ 

me  s 


A A 

n,+  = n,+ 

S 


The  relations  analogous  to  (B4)  are 

f - -jny  j*i  (z+z  ') 

n."  = e ' I A(r))e  e drj  ; 0 < z < z 


00  / 

- . -o  * A <*  -jiiy  - jn, (z+z  ) 

I1"T  = e'jPX  (B(ri)-A)e  e dri 

s J 


n /+  -jSx 

FI  = e 7 


J Qn). 


- jT»y  - j (z+z  1 ) 


■av  f - -j^y  -j^?z 
n"'  = e'jBX  , D(ri)e  e ^ dr] 

S v 


• 3x  p - -tty  j^2z 
n'  = e JBX  J E(T])e  e ^ dil 


v,Ta*iflearaBKfssp  -s.  ? 


Application  of  continuity  of  the  transverse  to  z fields  at  z = 0 yields 
via  ( B 1 ) the  following  equations: 


r , q\ 2-.  * " j^-i  z 

1 1 - ) j Be 

L \k'  J 


ve  - (!) ) ° 


[‘  - (D  ]■ 


- - VkJ  ) E 


-JHjZ  « -JJljZ  * a 

ate  h,  e C - Bqe  B = - uue  eH.,E  - Sr|D 

o 1 o 2 


-JKjZ  - 


C - uua  K . e 
o 1 


* -jH^Z  « - '■I 

B + 2uup  K.e  A=  -0riE  + ouu  k_D  ’ . 

o 1 o 2 


Elimination  of  D and  E between  (B21)  and  (B22)  yields 


A 2 2 
+ 2ujpok1A  = «»U01_k1(c-(-|)  ) + *2(l  - [£}  )]  B + 9ri(e-l)C 

2 2 _ a 

0 = - 3„(e-im  + »€0[k^e  -(■£))  + «2(l  - \|)  ) J C 


Similarly  to  (Bll)  the  system  determinant  A is 


2 

A = + k2(e  - (|)  ) (kj  + k2)  (h2  + €Hj) 


B = 2k^A 


».(«-  (f)V  - (if) 

X - (!)  X1  - (I)  X*1 + K2y(K2 + eKi) 


<-(! 


2TTU*Ue  - (!)  X1  - (!) 


_ r \ \k/  y-Vk, 

^ L k2  + £H1 


\k/  n 

+ — - — i 

h2  + K1 


t**>‘  J*#  * •*»  '*#*-•  * **  * ' 


' .•  • * • ■*  ’» 
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Appendix  C 

Determination  of  the  sheets  of  and  x.^  on  which  the  Sommerfeld 
poles  are  located. 


by 


The  poles  of  the  integrand  in  (14)  and  (15)  of  G and  H are  given 


^2+eXj  - 0 


Without  loss  of  generality  one  may  assume  that  the  poles  are 
located  on  the  proper  sheet  I x.^  < 0 of  x^  (otherwise  we  may  change 
signs  of  both  square  roots.  ) The  question  arises,  are  these  poles 
located  on  the  proper  sheet  of  Xj(I^  < 0 ) or  the  improper  sheet  of 
Xj  . Equation  C(l)  may  be  written: 


7k2-* 


e-a  + e 


e7k2-a2  = o , a2  = r]2+02 


C(l) 


C(2) 


H 

?! 

♦ 


4 

& 

Vf 


ft 


f 

w 

► 


and  are  given  by 


a2  = k2if 

p l + € 


Keeping  in  mind  that  ^ 0 , < 0 , we  have 


*2  = 'kJmr  ~ k/r 


so  that  on  the  proper  sheet  of  h ^ (x^  < ®)-  K2r  > ® 
In  a similar  fashion,  at  the  poles 


f 2 k2e 


Hi = Vk  ' irr = kv  i+t 


Since  with  e^>0  arri  < 0 , is  in  the  first  quadrant,  either  in 

the  first  or  the  third  quadrant  of  x.^.  The  former  corresponds  to  the 
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C(3) 


C(4) 


rt-f  ■*  **  v. 

- - ■ 


■ 'V'-^aSKrSSSF’i?'  * % 


T\ 


1 1'1 

V. 


improper  sheet  x^  > 0 and  since  in  the  first  quadrant,  x^  > 0 , The 

latter  yields  x. . < 0 >x.  <0.  To  determine  which  cases  are  possible 

’ li  lr 

let  us  rewrite  C(l)  in  terms  of  its  real  and  imaginary  parts. 


*2r+j  *2i  + <Cr+j€i)  (^r+jKli)  = ° 


C(5) 


or 


x,  +e  k,  -e. H,.  = 0 

2r  r lr  i li 


x0.  + e x..  + e.  x.  = 0 

2i  r li  l lr 


simultaneously 


C(6) 

C(7) 


The  following  signs  arise  in  C(6)  on  the  improper  sheet  of  x^ 


^2r+er  Klr-Gi  *11 


C(8) 


>0  >0  >0  <0  >0 


Thus,  all  terms  are  of  the  same  sign  and  cannot  add  up  to  zero  on  the 
improper  sheet  of  x^  . 

On  the  proper  sheet  of  X^  one  has  for  the  LFiS  of  C(6) 


x_  + € x - e . x. . 

2r  r lr  i ii 


>0  >0  < 0 <0  <0 


C(9) 


Thus,  the  first  term  is  positive  and  the  last  two  negative  and  balance  to 
zero.  To  verify  vanishing  of  LHS  of  C(7)  on  the  proper  sheet  of  x^  one 
has 

x.,.  + e x,.  + e.  x. 


2i 

<0 


r 'li  ' 1 lr 
>0  <0  <0  <0 


C(10) 


Thus,  the  first  two  terms  are  negative  and  the  third  positive.  The  above 
argument  establishes  the  location  of  Sommerfeld  poles  simultaneously  on 
the  proper  sheets  of  anc*  • 
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